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The Casimir force can be understood as resulting from the radiation pressure exerted by the 
vacuum fluctuations reflected by boundaries. We extend this local formulation to the case of par- 
tially transmitting boundaries by introducing reflectivity and transmittivity coefficients obeying 
conditions of unitarity, causality and high frequency transparency. We show that the divergences 
associated with the infiniteness of the vacuum energy do not appear in this approach. We give 
explicit expressions for the Casimir force which hold for any frequency dependent scattering and 
any temperature. The corresponding expressions for the Casimir energy are interpreted in terms of 
phase shifts. The known results are recovered at the limit of a perfect reflection. 
PACS: 03.65 - 42.50 - 12.20 



INTRODUCTION 



\q ' As a consequence of Heisenberg inequalities, the electromagnetic field exhibits quantum fluctuations even in the 
vacuum state. These "vacuum fluctuations" have been extensively studied in the domain of quantum optics these last 
years [Q. Casimir Q observed that the resulting vacuum energy {\fwj per field mode) depends upon the boundary 
conditions and particularly upon the positions of reflecting bodies. As a consequence, the vacuum fluctuations manifest 
themselves through macroscopic forces. 

These Casimir forces are usually computed by comparing the mode density in the absence and in the presence of 
perfectly reflecting boundaries H . In a local formulation, they can also be understood as resulting from the radiation 
pressure exerted by the vacuum fluctuations reflected by the boundaries . In the simplest case of a scalar field in 
a two-dimensional (2D) spacetime, one gets the following force between two pointlike mirrors at a distance q 

F=^ (1) 

Q., 

For the electromagnetic field in a four-dimensional (4D) spacetime, one gets the Casimir pressure (force measured per 
unit area) between two parallel plane and infinite mirrors || 

F - (0\ 



From now on, we use natural units where c = £o = fee — 1; however, we keep fi as a scale for vacuum fluctuations. 

A problem in any calculation of vacuum induced effects is to dispose of the divergences associated with the infinite- 
ness of the total vacuum energy. Usually, this is done by arbitrarily cutting off the mode density. The divergences of 
the energy in the absence and in the presence of the boundaries cancel each other, which leads to a finite net result at 
the limit of infinite cutoff frequency. In the local formulation, the field correlation function diverges when the fields 
are evaluated at the same point. A regularized stress tensor is obtained by splitting the two points and ignoring the 
divergent part 

Clearly, a more natural regularization should be provided by studying partially transmitting mirrors. Any real 
mirrors are certainly transparent at high frequencies so that the expression of the force should be regular. Using the 
physical properties of a dielectric constant, Lifshitz J8| has obtained regular expressions for the Casimir- Polder forces 
between macroscopic dielectric bodies P,[L0|. This approach is not free from problems, such as the effect of dissipation 
and the expression of energy density in a dispersive dielectric medium, which have been discussed in detail []TT ,p~2| . 

In this paper, we use a scattering approach to avoid these problems. We use the local formulation and intro- 
duce frequency dependent reflectivity and transmittivity coefficients for the two mirrors, which are assumed to obey 
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conditions of reality, unitarity, causality and high frequency transparency. The scattering coefficients determine the 
vacuum stress tensor and therefore the Casimir force on the mirrors. Explicit expressions, which come out as finite 
integrals, are obtained for any frequency dependent scattering. Eventually, the reflectivity function appears as a 
physical regulator and the known results (Q) and (|J) are recovered when the mirrors are totally reflecting over a large 
enough frequency interval. 

These results do not rely upon a detailed microscopic analysis of the scattering process and are not limited to 
dielectric mirrors. Dissipative processes are disregarded (unitary scattering) and all field fluctuations originate from 
the input vacuum. We also give the results for a non zero temperature input state. The Casimir energy can be deduced 
by integrating the force and is different from the integrated field energy when the reflectivities are frequency dependent. 
This expression, interpreted in terms of phase shifts, makes the connection with the mode density approach. 

A first part details the simple case where two pointlike mirrors are placed in the vacuum state, or in a thermal 
state, of a scalar field in a 2D spacetime. We come then to the problem of parallel plane and infinite mirrors which 
scatter the electromagnetic field in a 4D spacetime. Following Lifshitz Jp|, the effect of evanescent waves is accounted 
for. Expressions of the forces as integrals over imaginary frequencies || are given in an appendix. 

The following convention is used throughout the paper: any function f(t) defined in the time domain and its Fourier 
transform f[tu] are related through [] 
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ONE MIRROR IN 2D VACUUM 

In a 2D spacetime (one time coordinate t, one space coordinate x) , a free field is the sum of two counterpropagating 
fields 

$ = ip (t - x) + rjj (t + x) 
It can also be described by a twofold column matrix in the frequency domain 



$[w] = 



i/>[w] 



We consider now that the field is scattered by a mirror located at point q 

$ = 6{q - x) (<p- m (t - x) + tp out (t + x))+6{x- q) (t/j in (t + x)+ (p out (t - x)) 
where ip- ln and ip- m are the input fields, cp on t and ip on t the output fields (see Figure 1). 



out 



Pout 

► 



FIG. 1. One mirror scatters the two counterpropagating fields. 
With the hypothesis of a perfect reflection, the field vanishes on the mirror 

¥>out(*) = —1pia(t + 2q) 1paat(t) = —fm(t - 2<j) 

Equivalently 

VoutM = -ipm[uj]e~ 2 ' lujq Vw[H = -(pin[uj}e 2lujq 



1 The notation used in the original paper for Fourier transforms has been changed to a more convenient one. 
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In the general case of a partially transmitting mirror, the scattering of the field is described by a frequency dependent 
S matrix 

-2iujq 



4\, 



SM# ta [w] S[LU] 



s[oj\ r\u)\e 



(3) 



We assume that the S matrix obeys the following conditions. The values of S are real in the temporal domain 

s[—u] = s[u}* r[-u] = r[uj\* (4) 

The scattering is causal |H| 

s(t) = r(t) = for t < (5) 
s[uj] and r[u>] are analytic and regular for Imw > 
The scattering matrix is unitary (dissipation inside the mirror is neglected) 

S[w]S[u]^ = 1 

\s[oj]\ 2 + \r[uj}\ 2 = 1 s[w]r[<J\* + r[co]s[u]* = (6) 

The mirror is transparent at high frequencies and the reflectivity decreases sufficiently rapidly 

s[uj] — ► 1 u> \r[io]\ — > for lu — > oo (7) 

Because of the transparency condition (^), the reflectivity will provide a natural regulator to the ultraviolet diver- 
gence of the vacuum energy. The perfectly reflecting mirror corresponds to r = — 1 and s = at all frequencies and 
does not obey this condition. So, it will be preferable to consider the perfect mirror as the limit of a model obeying 
the conditions This amounts to consider the effect of non zero reflection delays and to let them go to zero at 

the end of the calculations. 

Quantum field theory provides the field correlation functions corresponding to the vacuum state. We will use the 
symmetric correlation functions which are the mean values of the anticommutators 

c(i) 
2 

+ r), <£/(*)}) =0 

<p'(t) = dt<p(t) 



<V(* + r), ¥/(«)}) =(W(t + r)^\t)}) 



(8) 



iuiiui' ({ <£>[<*>], tp{u> )}) = iujiuj' ({^[lo], ip(u> )}) = 2irS(u> 
Kjiu/({</>[w],^(a/)}) = 



with 



c(r) = lim — 
c[w] = ft |w| 



f 



1 



(e + ir) (e — ir)" 



(9) 

As usual in the local formulation, one evaluates the force as the difference between the radiation pressures exerted 
upon the left and right sides of the mirror. In a 2D spacetime, the component T xx of the stress tensor is equal to the 
energy density and one gets 



F = e L - e R 

e L = (^ n 2 (< -q) + + l)) e R - <^' n 2 (< + q) + ip'J t (t - q)) 

Substituting the output fields in terms of the input fields and using the unitarity condition, one obtains 

<^ u 2 t(i-<z)) - + = (v'-Jit-q)) = (rt n 2 (t + q)) - ^ 

Hence 



(10) 



' L = ' l: = ' 2tt C[UJ 



(11) 



The radiation pressures exerted upon the two sides of the mirror are infinite. However, they cancel each other exactly 
and the net force is zero. This is connected to the property that the impulsion density is zero in the vacuum state 
of a free field while the energy density is infinite. Note that the cancellation between infinite quantities is obtained 
inside a single physical situation and not by comparing two different situations. 
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TWO MIRRORS IN 2D VACUUM 



We study now the situation where two pointlike mirrors are placed in the 2D vacuum at positions q\ and q 2 (see 
Fig.2). 



IP 



out 



<p 



out 



FIG. 2. Two mirrors scatter the two counterpropagating fields. 

As is usual in the theory of a Fabry Perot cavity, the global scattering may be deduced from the elementary 
scattering matrices given by (||) and associated with the two mirrors 



<AiavM 

Vw M 

^out M 
J cav [ 



8i[w] 



s 2 [u>\ r 2 [uj\e 

o 2lL0q 2 



-2iuq2 



">PciivM 

</?cavM 
V'inM 



One can solve these equations to express the outcoming and the intracavity fields in terms of the input ones 

$ cav H = R[u]$ in [uj] 
The global scattering matrix S[u>] and the resonance matrix R[u>] are given by 

Si[w]s 2 M 



Su[u>] = 522 [w] = 
Si2[w] = ra[w]e- fa "- 
S2iM = ri[w]e- iw «- 



d[u 



s 2 b] 2 nMe iu « 



d[u 



siM 2 r 2 Me iw « 



r/[u 



(12) 



R 21 [u>] = 



d[w] 



R 12 [u] = 



s2H r iH e 



2iujq2 



d[uj] 



d[u] 
R22 M = 



52 M 



where 



d[uj] = 1 — r[ijj]e 2lW9 r[w] = ri[w]r2[w] 



q = 92 - qi 



(13) 



The matrices 5 and i? are real in the time domain, they are causal functions and the matrix 5 is unitary. 

Substituting the output fields in terms of the input fields and using the unitarity of the matrix S[w], one deduces 
that the energy densities and eR, again defined by (^), are still given by (^). But the intracavity energy density 
e cav is different 



(<t? 2 + 1b 1 2 ) 

\rcav 1 V-^cav/ 



dw 
2^ 



5M 



I Jill 



\Ri 



\R21 1 



|i?2 



1 - |r[w] 



|1 - r[w]e 2iw «| 



The function (calculated from equation (|l3|)) describes an enhancement (resp. suppression) of vacuum fluctuations 
seen inside the cavity at frequencies inside (resp. outside) the Airy peaks |l4|,[l^]. One deduces the Casimir force as 
the difference between the outer and inner radiation pressures 
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Fi = e L - e cav = F F 2 = e cav - e R = -F 



(14) 



F= / — c [u](l-g[w]) 
Jo 2lT 



(15) 



In the standard derivation of the Casimir force, one considers mirrors which are perfectly reflecting at all frequencies. 
In this case, the Airy peaks reduce to Dirac distributions 

1 - flM = 1 - - y^(w - n-) 
q q 

n 

The integral ( p"5| ) is then regularized by introducing a cutoff in the mode density and computed by applying the 
Euler-Mac Laurin summation rule Q 

Ba 7T 4 „, \ UijJ 



j=0 



where the B 2 k are the Bernouilli numbers B 2 = \] B4 = — .... This leads to the expression (|]) of the Casimir 
force. 

However, the introduction of a cutoff in the mode density is not necessary since the integral (^) is finite as soon 
as the reflectivity functions obey the conditions (0-0). One can write 



b = \ —c w - r : „. h c.c. (16) 

As c[w] is an even function of w, one deduces from the reality condition (^) 

F=rpc[ U ] 2 " rN rl (17) 

F = J2 F {e) F (e) = ~ [°° ^e 2tt ^/HcH (18) 

where each term F^' corresponds to a number i of roundtrips inside the cavity. This term can be transformed into 
a convolution product in the time domain, more precisely into an integral over I reflection delays t\, ... ti 

poo poo 

= - / dhrih) . . . / dt t r(U)c (2£q + h + ... + t e ) (19) 



One obtains the force as a sum 



Jo 



The reflection delays are positive and the two point correlation function c is evaluated at distances between the 
points greater than 2q. The expression of the force comes out free from the divergences usually associated with the 
infiniteness of the vacuum energy and the regulator e may be omitted in the expression (^) of c. 

In the limiting case where the reflection delays (ti, ... t£) are much shorter than the distance q between the two 
mirrors (large distance approximation), one can neglect them in the function c(2£q + 1\ + . . . + tg). It follows that the 
force depends only upon the reflectivities evaluated at zero frequency 



F = -J2 r o< 2e ^ r = r[0]=n[0]ra[0] 



i>0 

that is 



F = 4^(2) ( 2 °) 



C(p) = E|t ( 21 ) 
1=1 
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If the two mirrors are perfectly reflecting at zero frequency, one finds the expression ([j]) for the Casimir force 

(ci(2) = £). 

The equation ( p"8| ) allows one to evaluate the correction to the limit of ideal mirrors. It can be written as a formal 
series 

= A(_i( 2 ^- i ^)- 1 ( W rM £ )} +c.c. 

= A{( 2 ^-^)- 2 r[uf ) +c.c. (22) 

When r[w] is a very flat function, one recovers the large distance approximation where the force depends upon the 
value ro evaluated at zero frequency. The terms of the formal series which depend upon the derivatives of r[u] may 
be considered as corrections to this approximation. 

These discussions show that the reflectivity function r[u>] plays the same role as a regulator for the mode density. 
However, it is a regulator obeying the conditions (§-0) which cannot be real for real frequencies. We show in the 
appendix that it is possible to write the force as an integral over imaginary frequencies and that r[iu>] appears as a 
real regulator ||. 



CASIMIR ENERGY 



The local formulation led us to the expression of the force. One can deduce a Casimir energy U by integrating this 
force (see equation ([Til)) 

dU = -Fidgi - F 2 dq 2 = Fdq (23) 

We want now to compute this Casimir energy and to compare it with the integrated field energy Q . 
One deduces from equation (|l6| ) 

In ZTT 2 



(25) 



l-rwe 2,u « 
AM = «Log -V — ^ — 

One deduces respectively from the reality condition (^) and from the transparency condition (^) 

A[0] = wA[w] ^ for u -> oo 

so that the energy can be written 

u = f~^] (26) 

This is the standard form for the phase shift representation of the Casimir energy || . Actually, one checks from the 
expression ( |l2| ) of the S matrix that 

detS[w] = detS 1 [oj]detS 2 [Lu]e lAluj] 

where Sx[u>] and S 2 [uj] are the elementary S matrices associated with the two mirrors; A[w] is the q— dependent part 
of the sum of the two phase shifts at frequency u). 
Now, one gets from equation (|2~5|) 



{2qr[u] -id„r[w])e 2i »« 

O u A w = , : „■ h C.C. 

1 - r[uj}e 2luJ i 

= - (1 - g[u]) (2q + 8 U S) - g[uj] sin {2uoq + S) c^Log (l - p 2 ) (27) 

where S and p are the frequency dependent phase and modulus of the reflectivity function r. Usually, the Casimir 
effect is computed for frequency independent reflectivities. In this case, one deduces that the Casimir energy U is 
identical to the integrated field energy Q 

U = q (e cav - e L ) = -qF 

This result is compatible with the definition ( p3| ) of U since the force scales as q~ 2 in this case 0. But the reflectivities 
have to be frequency dependent in order to fulfill the transparency condition (^) . The resulting terms in the Casimir 
energy can be considered as corrections associated with the reflection delays upon the mirrors. 
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CASIMIR FORCE FOR NON ZERO TEMPERATURES 



Now, we derive the expression of the Casimir force for non zero temperatures MM ■ We suppose that the reflectivities 
are independent of the temperature. Therefore, we have only to change the correlation function of the input fields. 
The Casimir force is given by the integrals ( |i"5| ) or ( |l7| ) with a modified function c 

cr[u>] = Co[oj]n T [u] n T [uj] = 7TTT\ ( 28 ) 

tanh i -ffi J 

where cq\lu] corresponds to the zero temperature and nr[w] is the mean number of thermal photons per mode. The 
equations (19) are still valid with a function Cy[cj] which is the Fourier transform of 

..ha 2 irT 

Ct(T) = 2 a — ~ T~ 

7r sinh (ar) n 

It is known that the correlation function ct is the same as for an accelerated vacuum which also suggests 

that there are relations between the Casimir effect for non zero temperatures and the physics of black holes jl9| . 
At the large distance approximation, one gets 

F = -5>£c T (2^) 

e>o 

At the low temperature limit, one has the same results as previously. At the high temperature limit, one obtains an 
exponentially small force since 

Ct(t) „ -^T 2 exp (-^) for Tr » h 



T 1 ex v{-—) 

The thermal contribution exactly cancels the vacuum contribution for temperatures T higher than - . 

In the general case (no approximation on the distance), we obtain the Casimir free energy T by integrating the 
force j| 

F = dgHq ,T) F=r^^n T [u] (29) 



where A[u>] is still given by (|25h. Noting that (S is the entropy) 

U = T + TS = T - Td T T 
iodujiixluj] = —TOttix^} 



one also obtains 



An integration by parts then leads to the modified phase shift representation 

w = /"£v itM ^ am (30) 



Using equation fl27|), one reaches the same conclusions as previously about the comparison between the Casimir energy 
and the integrated field energy. 



TWO MIRRORS IN THE ELECTROMAGNETIC VACUUM 



We consider now that two parallel, plane and infinite mirrors are placed in the vacuum state of the electromagnetic 
field (plate separation q along the x direction). The scattering on each mirror is still described by a S matrix which 
obeys conditions of reality, unitarity, causality and high frequency transparency. Now, the reflection coefficients 
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associated with the two mirrors depend not only upon the frequency, but also upon the propagation direction and 
the polarization. We will denote 



7 = COS if K — LO COS if — LUJ K — LO sin if 

where tp is the incidence angle, n the wavevector along the x direction and k' the length of the transverse wavevector. 

The calculations of the Casimir force (measured per unit area) are mostly the same as in the 2D case, but one 
obtains it as a sum over the two polarizations p\ and p2 

F — F + F 

Each contribution F p can be written as an integral over the modes 

Fp = J J ^ 2 £ M c M(l-5 P k^])7 2 (31) 

The function c is the same as in the 2D case. The function g p is computed as previously, but its resonances are now 
determined by the wavevector measured along x 

r I l-Mw,K]| 2 

g p \uj,K\ = * 

\l~r p [Lu, K ]e 2 ^\ 2 

t p [lo,k] = r lp [u), K\r 2p [uJ, k] 

where r\ p and V2 P are the reflectivities associated with the two mirrors. In comparison with the 2D computation, 
there is in equation ( pl| ) an extra geometrical factor j 2 which is the ratio between the component T xx of the stress 
tensor (which gives the radiation pressure upon the mirror) and the energy density |20| ] . 

If the mirrors are supposed perfectly reflecting at all frequencies, Dirac distributions appear in the formula ( pi] ) 
which can be computed by a regularization and the application of the Euler-Mac Laurin summation rule Q 

The constant K goes to infinity with the cutoff frequency. But the even derivatives do not appear in the Euler-Mac 
Laurin sum and this constant is eliminated, which leads to expression (|^) of the force. 
In the general case of a frequency dependent reflectivity, one writes 

F * = J J 2^ HCM " 7 e-^-r p [^ 1 ] +C - C - (33) 

Now, one must specify more precisely the integration path for 7. Without further investigation, it could be admitted 
that it coincides with the segment [0, 1[. Using the same methods as for the 2D computation, one shows however 
that this integration path leads to a Casimir force which diverges at the limit of a perfectly reflecting mirror: there 
remains a contribution to the force analogous to the constant K appearing in (p^). 

Actually, the integration path for 7 must also include the upper part of the imaginary axis ]+ioo,0]. Imaginary 
values of 7 not too far from correspond to values of k' greater than, but not too much greater than u>. They are 
associated with vacuum fluctuations which enter the system through the side boundaries of the mirror and which 
contribute to the field outside the mirrors as evanescent waves. We will suppose that r p [ui, k] may be analytically 
continued plj for values of k in the quadrant [ReK > 0,lmft > 0] and also that the contribution of the evanescent 
waves to the force is given by the analytical continuation of the ordinary expression P2]| . Imaginary values of 7 
farther from correspond to larger values of k'. They are associated to "closed" input ports and have actually a zero 
contribution to the force. Assuming that |r p [o;, k] \ decreases sufficiently fast for large values of k and that |23j] 

\r p [u>, k]\ < 1 for to real and [ReK > 0,lmc > 0] 

one moves the integration path for 7 to ]1, 00 [ so that 



DC: 



duj [°° d7^ r A 2 2 r p [uj,uj-f} 



F p = — ricMtdV 9- ' T + cc. (34) 

This result is regular at the limit of perfect reflection. As for the 2D computation, one writes the force ( |34| ) as a 
sum over a number £ of roundtrips 
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t>o J ^ Z7r 



4%] = jf^k^ (35) 
(we have used a reality condition r p [o;, ft]* = r p [—L0 1 —k]). This expression can be written as a formal series 

h 



At the limit of a large separation, we retain the terms which depend upon the reflectivities evaluated around zero 
frequency but not those which depend upon the derivatives of the reflectivities. As 



we write 



R p ] ^ ~ ^ r o r o = {r p [n, k}} 

The value ro is evaluated at zero frequency and at normal incidence. This value is the same for the two polarizations 
and one obtains finally the total Casimir force as 

F = Y,r e C(2£q) C(r) = (36) 

that is 

^= 8^*6,(4) (37) 

where Cr(p) is defined by (^l|). At the limit of perfectly reflecting mirrors, one recovers the expression (||) for the 
Casimir force (Ci(4) = fo)- 

As in the 2D case, the expressions of the force can be written in terms of the reflectivities evaluated for imaginary 
frequencies (see the appendix). 

A phase shift representation can be obtained for any distance by integrating the force. One writes the energy U as 
a sum over the ordinary and evanescent waves 

JO 2lT J too 2lT 2 

where A [a;, k] is obtained by summing all the phase shifts at frequency u> and wavevector k 

A[w, k] = A P1 [u, k] + A P2 [uj, k] 

Moving the integration path leads to expressions in terms of the unphysical waves 
These expressions can be written in terms of 9 w A[w,a;7] by an integration by parts 

r°° duj r 1 d 7 hjLu 

l 2^J loo 2^~ 



/■°° dw r 00 d 7 ^ 7 w 3 . OAT 

Jo 2vr yj 2vr 6 
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Writing 9 W A[u>, ljj] as in equation (p7|), one obtains the Casimir energy as the sum of the integrated field energy and 
of corrections associated with the reflection delays upon the mirrors. The integrated field energy is here 



qF 



In the limit of a perfect reflection, this relation is compatible with the definition (|2^) of U [0 since the force scales as 



T 4 



The force for non zero temperatures is given by equations (^3j-|3J) with the function c modified according to equation 
(p8|). It is also given by equation (p5|) with 



At the large distance approximation, one obtains 



F = J2r e C(2lq) C(r) = -d 2 T 



r i« r i 
— r p [u;,K\ n T [u;\ 

ITT 



h „, a nT 



a 



7r 2 T tanh(ar) Ti 

The low temperature limit gives the known result. For high temperatures, one obtains by omitting exponentially 
small terms 



2T 

C(t) « — for Tt > h 

ITT 6 

F =4^sCr (3) (38) 



where Cx{p) is defined by (21). High temperatures correspond to a classical limit so that the expression does not 
contain h. For a perfect mirror, one recovers the known result Q. 

In the general case, one also obtains a phase shift representation of the free energy 

JO 27T Jtoo -™ 2 

duo f°° dj H^lu 2 
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APPENDIX A: EXPRESSIONS OF THE FORCES AS INTEGRALS OVER IMAGINARY FREQUENCIES 

Following the computation by Lifshitz of the Casimir-Polder forces ||, one can write the forces as integrals over 
imaginary frequencies. 

In the 2D case, one first shows that \r[uS\\ remains smaller than 1 for any to in the upper half plane Imw > 0. 
The modulus \r[uj] \ is less than 1 on the real axis (from the unitarity and reality conditions); it goes to for large 
values of oj (from the transparency condition); as r[oj] is analytic for lmuj > 0, the conclusion is a consequence of the 
Phragmen-Lindelof theorem |l3| . 

Therefore, the integrand appearing in (^) has no poles in the upper half plane. Using the transparency condition, 
one can then move the integration path [0, oo[ to the upper part of the imaginary axis [0,ioo[ 

huj r[iuj] 
dui- 7, f — t + c.c. 

2ir e 2u i - r[iuj] 

Furthermore, r[iuj] is real since it is the Laplace transform of the causal and real function r 

dt r(i)e" wt 
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It follows that the Casimir force is 



F = 



dto 



^2uq 



This equation can conveniently be used for evaluating the Casimir force for any reflectivity functions 



F = 



4irq 2 



du- 



ur 


iu 








e u _ 


r 


iu 






2 1. 



(Al) 



At the large separation limit, r ^ may be replaced by ro and one obtains an integral expression equal to the series 
of equation (|20|) 



F = 



The integral (Al) can be written as a formal series 



du- 



ur 



47r<z Jo e u — tq 



F = J2F (e) 

F W = ll(2eq-d u )- 1 wr[iu] l \ = - \(2£q - c^)" 2 r[i<j] e \ 

This equation and the expression ( |22| ) are formally equivalent; however, they have a different status: r[iu>] is real 
and plays the role of a regulator for the mode density while r[u>], which appears in (|22]), cannot be real because of 
causality and high frequency transparency. 

One can also write the 4D Casimir force as an integral over imaginary frequencies. We start from the integral 
( |33| ) where the integration path for 7 includes the ordinary contributions (0 < 7 < 1) as well as the evanescent 
waves (7 imaginary) and we treat separately these two contributions. For the ordinary contributions, it follows from 
causality that: r p [w,w7] is analytic in u> for [Reuj > 0,lmcij > 0] and 7 a fixed real number; |r p [o;, W7]| < 1 for 
[Rew > 0,lmw > 0] and 7 real. Using the high frequency transparency, we move the integration path for u> to the 
upper imaginary axis (with 7 kept constant). For the evanescent waves, we use the analogous properties j22j: t p [lu, k] 
is analytic in lu for [Reuj > 0,lmw > 0] and k a fixed imaginary number; |r p [w,K]| < 1 for [Reuj > 0,lmcj > 0] and 
k imaginary. Using the high frequency transparency, we move the integration path for u> to the upper imaginary axis 
(with k kept constant). Finally, we use the property: r p [u, k] real for ui and n imaginary. 

At the end of the process, one obtains 



p — 



dn k 2 



dto- 



r„\iu>, ik\ 



,2nq 



r n nu}, ik\ 



At the large separation limit, r v \iuj, in] is replaced by a constant r$ and one obtains the Lifshitz' integral expression 



equal to the series of equation ( |37| ) 



F 



du- 



u 3 r 



'0 
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